arXiv: 1503.03197vl [astro-ph.SR] 11 Mar 2015 


Draft version March 12, 2015 

Preprint typeset using DTp;X style emulateapj v. 5/2/11 


OGLE-2013-BLG-0578L: MICROLENSING BINARY COMPOSED OF A BROWN DWARF AND AN M DWARF 

H. Park 1 , A. Udalski 2,5 , C. Han 1,6 , 

AND 

R. Poleski 2,3 , J. Skowron 2 , S. Kozlowski 2 , L. Wyrzykowski 2,4 , M. K. Szymanski 2 , P. Pietrukowicz 2 , G. Pietrzynski 2 , 

I. SOSZYNSKI 2 , K. ULACZYK 2 
(The OGLE Collaboration) 

'Department of Physics, Institute for Astrophysics, Chungbuk National University, Cheongju 371-763, Korea 
2 Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warszawa, Poland 
department of Astronomy, Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA and 
institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK 
Draft version March 12, 2015 

ABSTRACT 

Determining physical parameters of binary microlenses is hampered by the lack of information about the an¬ 
gular Einstein radius due to the difficulty of resolving caustic crossings. In this paper, we present the analysis 
of the binary microlensing event OGLE-2013-BLG-0578, for which the caustic exit was precisely predicted in 
advance from real-time analysis, enabling to densely resolve the caustic crossing and to measure the Einstein 
radius. From the mass measurement of the lens system based on the Einstein radius combined with the addi¬ 
tional information about the lens parallax, we identify that the lens is a binary that is composed of a late-type 
M-dwarf primary and a substellar brown-dwarf companion. The event demonstrates the capability of current 
real-time microlensing modeling and the usefulness of microlensing in detecting and characterizing faint or 
dark objects in the Galaxy. 

Subject headings: binaries: general, brown dwarfs - gravitational lensing: micro 


1. INTRODUCTION 

It is known that low-mass stars comprise a significant frac¬ 
tion of stars in the Solar neighborhood and the Galaxy as a 
whole. The Galaxy may be teeming with even smaller mass 
brown dwarfs. Therefore, studying the abundance and prop¬ 
erties of low-mass stars and brown dwarfs is of fundamen¬ 
tal importance. There hav e been surveys searching for very 
low-m ass (VLM) objects dReid et all 120081 lAberasturi et all 
120141) . but these surveys are limited to the immediate solar 
neighborhood. As a result, the sample of VLM objects is 
small despite their intrinsic numerosity and thus our under¬ 
standing about VLM objects is poor. 

Microlensing surveys detect objects through their gravita¬ 
tional fields rather than their radiation and thus microlensing 
can provide a powerful probe of VLM objects. However, the 
weakness of microlensing is that it is difficult to determine 
the lens mass for general microlensing events. This difficulty 
arises due to the fact that the time scale of an event, which is 
the only observable related to the physical lens parameters, re¬ 
sults from the combination of the lens mass, distance and the 
relative lens-source transverse speed. As a result, it is difficult 
to identify and characterize VLM objects although a signifi¬ 
cant fraction of lensing events are believed to be produced by 
these objects. 

However, it is possible to uniquely determine the physical 
lens parameters and thus identify VLM objects for a subset of 
lensing events produced by lenses composed of two masses. 
For unique determinations of the physical lens parameters, it 
is required to simultaneously measure the angular Einstein ra¬ 
dius 9e and the microlens parallax ite that are related to the 
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lens mass M and distance to the lens l)\ by 


Mm 


K7T E ’ 


d l = 


AU 

tteOe + its ’ 


(1) 


respectively (IGouldl l2000h . Here k = 4G/(c 2 AU), tt s = 
AU/Ds is the parallax of the source star, and D$ is the dis¬ 
tance to the source. The angular Einstein radius is esti¬ 
mated by analyzing deviations in lensing light curves caused 
by the finite size of the lensed source stars (IGouldl 119941: 
iNemiroff & Wickramasing he 1994: 1 Witt & Maoll 1 9941) . For a 
single-lens events, finite-source effects can be detected only 
for a very small fraction of extremely high-magnification 
events where the lens-source separation at the peak magni¬ 
fication is equivalent to the source size. On the other hand, 
light curves of binary-lens events usually result from caustic 
crossings during which finite-source effects become impor¬ 
tant and thus the chance to detect finite-source effects and 
measuring the Einstein radius is high. Furthermore, binary- 
lens events tend to have longer time scales than single-lens 
events and this also contributes to the higher chance to mea¬ 
sure the lens parallax. In fact, most known VLM lensing 
objects were identified through the channel of binary-lens 
events dH wang et alJl2010t lShin et alJl2012t lChoi et al.ll2013b 

lHan et alJl2013b iPark et alJl2013HJung et al1l2015l). 

Despite the usefulness of binary-lens events, the chance to 
identify VLM objects has been low. One important reason 
for the low chance is that caustic crossings last for very short 
periods of time. The duration of a caustic crossing is 
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where fg is the Einstein time scale of the event, <f> is the en¬ 
trance angle of the source star with respect to the caustic line, 
and p = 9*1 Be is the angular source radius 0* normalized to 
the angular Einstein radius 9e . Considering that the Einstein 
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time scale is ~ (0)10 days and the Einstein radius is ( O ) 
milli-arcsec for typical Galactic microlensing events, the du¬ 
ration of a caustic crossing is ~ (O) hours for Galactic bulge 
source stars with angular stellar radii of ~ ( O ) 1-10 /i-arcsec. 
Therefore, it is difficult to densely resolve caustic crossings 
from surveys that are being carried out with over hourly ob¬ 
servational cadences. 

Another reason for the low chance of resolving caustic 
crossings is the difficulty in predicting their occurrence. Caus¬ 
tics produced by a binary lens form a single or multiple sets 
of closed curves and thus caustic crossings always come in 
pairs. Although it is difficult to predict the first crossing (caus¬ 
tic entrance) based on the fraction of the light curve before the 
caustic entrance, the second crossing (caustic exit) is guaran¬ 
teed after the caustic entrance. To resolve the short-lasting 
caustic exit, it is required to precisely predict the time of the 
caustic crossing so that observation can be focused to resolve 
caustic crossings. This requires vigilant modeling of a lens- 
ing event conducted with the progress of the event followed 
by intensive follow-up observation. 

In this paper, we report the discovery of a VLM binary 
that was detected from the caustic-crossing binary-lens mi¬ 
crolensing event OGLE-2013-BLG-0578. The caustic exit of 
the event was precisely predicted by real-time modeling, en¬ 
abling dense resolution and complete coverage of the caustic 
crossing. Combined with the Einstein radius measured from 
the caustic-crossing part of the light curve and the lens par¬ 
allax measured from the long-term deviation induced by the 
Earth’s parallactic motion, we uniquely measure the lens mass 
and identify that the lens is a VLM binary composed of a low- 
mass star and a brown dwarf. 

2. OBSERVATION 

The microlensing event OGLE-2012-BLG-0578 occurred 
on a star located toward the Galactic Bulge direction. The 
equatorial coordinates of the lensed star are (a,5)j2ooo = 
(17 h 59 m 59 s .85, -29°44'06".9), that correspond to the Galactic 
coordinates (l,b) = (0°.90,-3°.10). The event was first noticed 
on April 22, 2013 from survey observations conducted by the 
Optic al Gravitational Lensing Experiment (OGLE: [Udalskil 
120031 ) using the 1.3m Warsaw telescope at the Las Campanas 
Observatory in Chile. Images were taken using primarily in 
/-band filter and some V-band images were also taken to con¬ 
strain the lensed star (source). 

In Figure Q] we present the light curve of the event. The 
light curve shows two distinctive spikes that are characteristic 
features of a caustic-crossing binary-lens event. The caus¬ 
tic crossings occurred at HJD' = HJD — 2450000 ~ 6426.0 
and 6461.8. Although the event was first noticed before 
the caustic crossings, the caustic entrance was missed. With 
the progress of the event, it became clear that the event was 
produced by a binary lens from the characteristic “U”-shape 
trough between the caustic crossings. 

Although the first caustic crossing was missed, the second 
crossing was densely resolved. Resolving the crossing be¬ 
came possible with the prediction of the caustic crossing from 
vigilant modeling of the light curve followed by intensive 
follow-up observation. It is known that reliable prediction 
of the second caustic crossing is difficult based on the light 
curve before the minimum between the two caustic crossings 
(iJarosz vnski & Mao 2001). With the emergence of the cor¬ 
rect model after passing the caustic trough, we focused on the 
prediction of the exact caustic-crossing time. The first alert 
of the caustic exit was issued on June 17, 2013, 1.3 days be- 
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Fig. 1.— Light curve of OGLE-2013-BLG-0578. Also plotted is the best- 
fit model (“orbit+parallax” model with mo >0). The two bottom panels show 
the residuals from the modeling with and without considering lens-orbital and 
parallax effects. The upper panel shows the enlargement of the caustic-exit 
region and the times of caustic-crossing alerts. 

fore the actual caustic exit. On the next day, the second alert 
was issued to predict more refined time of the caustic exit. In 
response to the alert, the OGLE experiment, which is usually 
operated in survey mode, entered into “following-up” mode 
observation by increasing the observation cadence. Thanks to 
the intensive follow-up observation, the caustic exit was com¬ 
pletely and densely covered. See the upper panel of FigureQ] 

In order to securely measure the baseline magnitude and de¬ 
tect possible higher-order effects, observation was continued 
after the caustic exit to the end of the Bulge season. From 
these observations, we obtain 2284 and 27 images taken in / 
and V bands, respectively. Photometry of the event was done 
by using the customized pipeline ( IUdalskifl2003l) t hat is based 
on the Difference Ima ging Analysis method ifAlard & Luptonl 
ll998tlWozniakl f2000). We note that the /-band data are used 
for light curve analysis, while the V-band data are used for the 
investigation of the source type. 

It is known that photometric errors estimated by an auto¬ 
matic pipeline are often underestimated and thus errors should 
be readjusted. We readjust error bars by 


s' = k(e 2 + e 2 min ) l/2 . 


(3) 


Here e mm is a term used to make the cumulative distribution 
function of \ 2 as a function of lensing magnification becomes 
linear. This process is needed to ensure that the dispersion of 
data points is consistent with error bars of the source bright¬ 
ness. The other term k is a scaling factor used to make \ 2 P er 
degree of freedom (dof) becomes unity. 


3. ANALYSIS 

We analyze the event by searching for the set of lensing 
parameters (lensing solution) that best describe the observed 
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light curve. Basic description of a binary-lens event requires 
seven standard lensing parameters. Three of these parameters 
describe the source-lens approach including the time of the 
closest approach of the source to a reference position of the 
lens, to, the lens-source separation at f 0 in units of the Einstein 
radius, no, and the time required for the source to cross the 
Einstein radius, ?e (Einstein time scale). In our analysis, we 
set the center of the mass of the binary lens as a reference po¬ 
sition in the lens plane. Another two parameters describe the 
binary lens including the projected binary separation in units 
of the Einstein radius, s (normalized separation), and the mass 
ratio between the lens components, q. Due to the asymmetry 
of the gravitational held around the binary lens, it is needed 
to define the angle between the source trajectory and the bi¬ 
nary axis, a (source trajectory angle). The last parameter is 
the normalized source radius p, which is needed to describe 
the caustic-crossing parts of the light curve that are affected 
by finite-source effects. 

It is often needed to consider higher-order effects in order 
to precisely describe lensing light curves and this requires to 
include additional lensing parameters. For long time-scale 
events, such effects are caused by the positional change of an 
observer induced by the orbital motion of the Earth around the 
Sun (“parallax effect”: iGouldll 1992i ) and/or the change of the 
binary separation and orientatio n caused by the orbital motion 
of the lens (“lens orbital effect’ Y lDominikil 998tlAlbrow et aTl 
120001) . The analyzed event lasted throughout the whole Bulge 
season and thus these effects can be important. The paral¬ 
lax effect is described by two parameters, tte,n and ~e.e, that 
are the two components of the lens parallax vector 7 Te pro¬ 
jected onto the sky along the north and east equatorial coor¬ 
dinates, respectively. The direction of the lens parallax vec¬ 
tor corresponds to the relative lens-source proper motion and 
its magnitude corresponds to the relative lens-source parallax 
7 r re ] = ALY/Of" 1 -Dg 1 ) scaled to the Einstein radius of the lens, 

t>E 


To the first-order approximation, the lens orbital effect is 
described by two parameters, ds/dt and da/dt , which are 
the change rates of the normalized binary separation and the 
source trajectory angle, respectively. 

To model caustic-crossing parts of the light curve, it is 
needed to compute magnifications affected by finite-source 
effects. To compute finite-source magnifications, we use 
the numerical method of the inverse ray - shooti ng technique 
dKavser et al.l Il986t ISchneider & Weissl 119861) in the im¬ 
mediate neighboring region around caustics and the semi- 
analytic hexadecapole approximation dPeicha & Hevrovskvi 
2009l lGouidl l2008) in the outer region surrounding caustics. 
We consider the effects of the surface brightness variation of 
the source star. The surface brightness is modeled as 


Saoc1-T a 



(5) 


where T A is the linear limb-darkening coefficient, A is the 
passband, and ip is the angle between the line of sight to¬ 
ward the source star and the normal to the source surface. We 
adopt the limb-darkening coefficients (IV, [’)) = (0.62,0.45) 
from lClaretl (2000 ) based on the source type. The source type 
is determined based on its de-reddened color and brightness. 
See Section 4 for details about how the source type is deter¬ 
mined. 



Fig. 2.— The source trajectory (the curve with an arrow) with respect to 
the caustic (the closed curve with 6 cusps). The two small dots marked by 
M\ and M 2 represent the positions of the binary lens components. All lengths 
are scaled by the angular Einstein radius corresponding to the total mass of 
the binary lens. Due to the change of the relative positions of the binary 
lens components caused by the orbital motion, the caustic varies in time. We 
present 4 caustics and lens positions corresponding to the times marked in 
right upper comer of the panel. The source trajectory is curved due to the 
combination of the lens-orbital and parallax effects. 


Searching for the best-fit solution of the lensing parameters 
is carried out based on the combination of grid-search and 
downhill approaches. We set (s,q,a) as grid parameters be¬ 
cause lensing magnifications can vary dramatically with small 
changes of these parameters. On the other hand, magnifica¬ 
tions vary smoothly with changes of the remaining param¬ 
eters, and thus we search for these parameters by using a 
downhill approach. We use the Markov Chain Monte Carlo 
(MCMC) method for the downhill approach. Searching for 
solutions throughout the grid-parameter spaces is important 
because it enables one to check the possible existence of de¬ 
generate solutions where different combinations of the lensing 
parameters result in similar light curves. 

In the initial search for solutions, we conduct modeling of 
the light curve based on the 7 basic binary-lensing parameters 
(“standard model”). From this, it is found that the event was 
produced by a binary with a projected separation very close 
to the Einstein radius, i.e. s ~ 1.0. Caustics for such a reso¬ 
nant binary form a single big closed curve with 6 cusps. The 
two spikes of the light curve were produced by the source 
trajectory passing diagonally through the caustic. See Fig¬ 
ure [2] where we present the source trajectory with respect to 
the causticQ The estimated mass ratio between the binary 
components is q ~ 0.2-0.3. 

Although the standard model provides a fit that matches the 
overall pattern of the observed light curve, it is found that 
there exists some residual in the region around HJD' ~ 6410. 
See the bottom panel of Figure Q] We, therefore, check 
whether higher-order effects improve fit. We find that sep¬ 
arate consideration of the parallax effect (“parallax model”) 

7 We note that the source trajectory is curved due to the combination of 
the lens-orbital and parallax effects. We also note that the caustic varies in 
time due to the positional change of the binary-lens components caused by 
the orbital motion. 
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TABLE 1 

Lensing parameters 


Parameter 

uo > 0 

«0 < 0 

X 2 /dof 

2303.1/2300 

2328.1/2300 

to (HJD) 

2456440.13 ±0.11 

2456439.99 ± 0.08 

u 0 

0.109 ± 0.004 

-0.120 ± 0.003 

t E (days) 

72.11 ±0.84 

70.79 ±0.10 

s 

1.027 ± 0.004 

1.035 ±0.002 

q 

0.260 ± 0.004 

0.241 ± 0.003 

a (rad) 

0.676 ± 0.007 

-0.647 ± 0.007 

P GO 3 ) 

0.96 ± 0.01 

0.97 ± 0.01 

™E,N 

0.54 ± 0.08 

-0.51 ±0.05 

7TE ,E 

-0.53 ± 0.03 

-0.72 ± 0.02 

ds/dt (yr” 1 ) 

0.48 ±0.14 

0.40 ± 0.03 

da/dt (rad yr -1 ) 

1.89 ±0.24 

-1.75 ±0.13 


and the lens-orbital effect (“orbital model”) improves fit by 
A\ 2 = 579 and 344, respectively. When both effects are si¬ 
multaneously considered (“orbit+parallax model”), the fit im¬ 
proves by Ay 2 = 616, implying that both effects are impor¬ 
tant. In Figure[3] we present the contours of Ay 2 in the space 
of the higher-order lensing parameters. Contours marked in 
different colors represent the regions with Ay; 2 < 1 (red), 4 
(yellow), 9 (green), 16 (sky blue), 25 (blue), and 36 (pur¬ 
ple). It shows that the higher-order effects are clearly de¬ 
tected. Considering the time gap between the caustic cross¬ 
ings that is approximately a month and the long duration of 
the event that lasted throughout the whole Bulge season, the 
importance of the higher-order effects is somewhat expected. 

It is known that lensing events with higher-order effects 
are subject to the degeneracy caused by the mirror sym¬ 
metry of the source trajectory with respect to the binary 
axis dSkowron et al.ll201ll) . This so-called “ecliptic degen¬ 
eracy” is important for Galactic Bulge events that occur 
near the ecliptic plane. The pair of the solutions resulting 
from this degeneracy have almost identical parameters except 
(uo,a,Tr^,da/dt) — > ~(uo,a,irE,N,da/dt). It is found that 
mo > 0 is marginally preferred over the mq < 0 solution by 
Ay; 2 = 25.0, which corresponds to formally ~ 5a level dif¬ 
ference. However, this level of Ay 2 can often occur due to 
systematics in data and thus one cannot completely rule out 
the mo < 0 solution. In Table Q] we present the best-fit pa¬ 
rameters of both mo > 0 and uq < 0 solutions. We note that the 
uncertainties of the lensing parameters are estimated based on 
the distributions of the parameters obtained from the MCMC 
chain of the solution. We also present the best-fit model light 
curve (mo > 0 solution) and the corresponding source trajec¬ 
tory with respect to the lens and caustic in Figures Q] and [2] 
respectively. 

4. PHYSICAL PARAMETERS 

By detecting both finite-source and parallax effects, one 
can measure the angular Einstein radius and the lens parallax, 
which are the two quantities needed to determine the mass 
and distance to the lens. The lens parallax is estimated by 
tie = (7r| n + n E fi) from the parallax parameters determined 
from light-curve modeling. 

In order to estimate the angular Einstein radius, it is needed 
to convert the measured normalized source radius p into 9e 
by using the angular radius of the source star, i.e. Oe = 
6*/p. The angular source radius is estimated based on the 
de-reddened color (V -1 )o and brightness Iq of the source 
star which are calibrated by using the centroid of the Bulge 



n “ ds/dt (yr->) da/dt (yr“‘) 


Fig. 3.— Contours of A\ 2 in the space of the parallax and lens-orbital pa¬ 
rameters for the best-fit model («o > 0 model). Contours marked in different 
colors represent the regions with A\ 2 < 1 (red), 4 (yellow), 9 (green), 16 
(sky blue), 25 (blue), and 36 (purple). 



Fig. 4.— Location of the lensed star with respect to the centroid of giant 
clump in the color-magnitude diagram of neighboring stars. 


giant clump on the color-magnitude diagram as a reference 
(lYoo et aP l2004). By adopting the color a nd brightness of the 
clump centroid (V -IJ)q c = (1.06,14.45) (iBensbv et al.l201 ll : 
iNataf et all 1 201 3l) . we estimate that (V-1,1) o = (0.72,17.68) 
for the source star, implying that the source star is a G-type 
subgiant. Figure |4] shows the locations of the source and the 
centroid of giant clump in the color-magnitude diagram of 
stars around the source star. We t hen translate V-I into V —K 
using the color-color relation of [Bessell & Brettj (119881) and 
obtain the angular source radius by using the relation between 
V — K and 0* of lKervella et al,l (120041) . The determined angu¬ 
lar source radius is 0* = 0.93 ± 0.06 pas. Then, the Einstein 
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TABLE 2 

Physical parameters 


of the obtained lensing solution, we compute the projected 
kinetic to potential energy ratio (KE/PE)_l by 


Parameter 

«o > 0 

Ho < 0 

Einstein radius (mas) 

0.97 ± 0.07 

0.96 ± 0.07 

Proper motion (mas yr -1 ) 

4.90 ± 0.35 

4.96 ± 0.35 

Total mass (Mq) 

0.156 ± 0.017 

0.133 ±0.011 

Mass of primary ( Mq) 

0.124 ±0.014 

0.107 ±0.009 

Mass of companion (Mq) 

0.032 ± 0.004 

0.026 ± 0.002 

Distance (kpc) 

1.16 ± 0.11 

1.02 ±0.08 

Projected separation (AU) 

1.16 ± 0.11 

1.02 ±0.08 

(KE/PE)j_ 

0.48 ± 0.20 

0.32 ± 0.07 


radius is estimated as 9e = 0.97 ± 0.07 mas for the best-fit 
solution (no > 0 model). We note that uq < 0 model results 
in a similar Einstein radius due to the similarity in the mea¬ 
sured values of p. Combined with the measured Einstein time 
scale Ie, the relative lens-source proper motion is estimated as 
p = 6e/ tE = 4.90 ± 0.35 mas yr -1 . 

With the measured lens parallax and Einstein radius, we de¬ 
termine the mass and the distance to the lens using Equation 
([]}• In Table [2] we list the physical parameters of the lens 
system corresponding to the uq > 0 and uq < 0 solutions. We 
note that the estimated parameters from the two solutions are 
similar. According to the estimated mass, the lens system is 
composed of a substellar brown-dwarf companion and a late- 
type M-dwarf primary. The distance to the lens is Z>l <1.2 
kpc and thus the lens is located in the Galactic disk. The pro¬ 
jected separation between the lens components is r± = sD\0e 
is slightly greater than 1 AU. In order to check the validity 


/KE\ _ (rj_/AU) 2 |7lM 2 (da\ 

(pe7_l “ 8tt 2 (M/Mo) V7<*7 V - -*/ 


( 6 ) 


(iDong et al.ll2009 ). The estimated ratio is (KE/PE)j_ < 1 for 
both mq > 0 and uq < 0 solutions and thus meets the condition 
of a bound system. 


5. CONCLUSION 

We presented the analysis of a caustic-crossing binary-lens 
microlensing event OGLE-2013-BLG-0578 that led to the 
discovery of a binary system composed of a substellar brown- 
dwarf companion and a late-type M-dwarf primary. Identifi¬ 
cation of the lens became possible due to the prediction of the 
caustic crossing from vigilant real-time modeling and resolu¬ 
tion of the caustic from prompt follow-up observation. The 
event demonstrates the capability of current real-time model¬ 
ing and the usefulness of microlensing in detecting and char¬ 
acterizing faint or dark objects in the Galaxy. 
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